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Abstract. We study the impact of a spatially homogeneous yet non-stationary 
dielectric permittivity on the dynamical and spectral properties of light. Our choice 
of potential is motivated by the interest in PT-symmetric systems as an extension of 
quantum mechanics. Because we consider a homogeneous and non-stationary medium, 
VT symmetry reduces to time-reversal symmetry in the presence of balanced gain 
and loss. We construct the instantaneous amplitude and angular frequency of waves 
within the framework of Maxwell’s equations and demonstrate the modulation of light 
amplification and attenuation associated with the well-defined temporal domains of 
gain and loss, respectively. Moreover, we predict the splitting of extrema of the angular 
frequency modulation and demonstrate the associated shrinkage of the modulation 
period. Our theory can be extended for investigating similar time-dependent effects 
with matter and acoustic waves in T^T-symmetric structures. 


Keywords: Amplitude and frequency modulations, VT symmetry in optics, time-varying 
wave propagation, analytical solutions to Maxwell’s equations 

1. Introduction 

During the past years, a new class of Hamiltonians has been widely investigated, which 
extends quantum mechanics from the Hermitian into the non-Hermitian (complex) 
domain [1]. Despite the lack of Hermiticity, Bender et al. have shown in their seminal 
papers that a Hamiltonian can have real eigenspectra if it possesses so-called parity¬ 
time {VT) symmetry [SIE]. Such a symmetry means there is invariance of the theory 
under parity (spatial) reflection V: p ^ —p, x —x, and time reflection T: p —)■ —p 
{t —t), i —i, X ^ X, where p and x are the momentum and position operators, 
respectively, while t is the time coordinate and i is the imaginary unit. This combined 
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VT symmetry leads to subtle changes in the unitary evolution of the system and 
modihcation of the inner product giEiiniiT]. As VT symmetry represents an extension of 
quantum mechanics, it is nowadays used in various different contexts, such as quantum 
reflection iiiniiio] and chaos m, and has even been generalized to fermionic [HI [13], 
gyrotropic [mils] and magnetic systems ng. 

Although the concept of VT symmetry was originally introduced in quantum 
mechanical systems, one has found experimental evidence and also a wide range of 
applications in classical optics. In 2010, Riiter et al. were the hrst to realize a PT-optical 
coupled system that involves well-dehned regions with gain and loss regimes, inherent to 
the complex-valued refractive index na. Such an extension of the concept of spacetime 
reflection into the classical domain stems from the works of El-Ganainy et al. [18] and 
Makris et al. [191EO] , who have employed the similarity between the Schrodinger and a 
scalar approximation of Maxwell’s equations to describe the dynamics of light beams in 
VT -symmetric optical lattices. There have been further theoretical [211 ESI [23l [23] and 
experimental [251 EHl ET] studies dealing with the implementation of the parity-time 
reversal symmetry in optics especially relevant for the development of new artihcial 
structures and materials [see also the recent review paper [2H] and references cited 
therein]. 

PT-symmetric structures have mainly been investigated in the spatial domain 
(that is, for time-independent complex potentials) and little attention has been paid 
to the study of the non-stationary regime. The signihcance of consideration of the 
time-dependent potentials both in the quantum [291 EQl EH [321 E3] and classical 
theories [SH [35] arises from the attempt to examine the full time evolution of the system. 
Despite the recent works, however, there are no rigorous analytical studies of wave 
equations with non-stationary complex potentials possessing time reflection symmetry. 
Given that the energy (frequency) and time are conjugate variables, a successful solution 
of PT-symmetric time-dependent Maxwell’s equations would constitute a complete 
characterization of the dynamical and spectral features of light. 

The purpose of this paper is to study the dynamics of light in time-dependent 
optical potentials having VT symmetry. We consider a spatially homogenenous system, 
for which VT reduces to symmetry under time reversal, albeit in the presence of both 
gain and loss. In view of this, we calculate both the instantaneous amplitude and 
angular frequency of waves and show how the complex-valued dielectric permittivity 
controls the light in the temporal domain. The resulting amplihcation and attenuation 
of the amplitude is demonstrated to be associated with well-dehned regimes of gain 
and loss, respectively. A comparison with the angular frequency modulation by a real 
permittivity is provided to reveal the impact of PT-symmetric potentials in that we 
observe (i) splitting of the extrema and (ii) a shrinkage of the frequency modulation 
period. Moreover, a particular emphasis is placed on studying modulations of light 
amplihcation and attenuation for experimentally accessible values of the modulated 
permittivity. 

The paper is organized as follows. In Section [2] we briehy discuss VT symmetry 
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in optics by considering space-independent but time-varying modulations of a complex¬ 
valued dielectric permittivity. For such a modulated permittivity, symmetric under 
time-reversal, we construct an analytical solution to Maxwell’s equations when the 
modulation rate is assumed to be much smaller than the wave frequency. This solution 
is further exploited in Section |3] to derive and analyze the instantaneous amplitude 
and angular frequency of light. Peculiar properties, such as modulations of amplitude 
amplihcation and attenuation as well as the splitting of frequency extrema and shrinkage 
of frequency modulation, are discussed in detail. Finally, conclusions are given in 
Section HI 

2. Solution to Maxwell’s equation for space-independent but time-varying 
dielectric permittivity obeying VT symmetry 

We start with a brief discussion of general properties of classical optical systems 
possessing VT symmetry. Building on the formal equivalence of the Schrodinger 
equation with the paraxial Helmholtz (Maxwell) equation, we identify the complex 
refractive index n = 3? (n)-l-iQ' (n) as the optical potential, the real (3? (n)) and imaginary 
(Q'(n)) parts of which are, correspondingly, even and odd functions of spacetime 
coordinates to ensure the VT invariance of the theory [m [ig [IS EU]. Likewise, since 
for non-magnetic structures the real (3? (e)) and imaginary (9= (e)) parts of the dielectric 
permittivity are dehned via 3? (e) = [3f?(n)]^ — [Q=(n)]^ and O'(e) = 23? (n) 3= (n), the 
symmetry, 3? [e (x, t)] = 3? [e (—x, —t)], and anti-symmetry, 9= [e (x, t)] = —O' [e (—x, —t)], 
relations guarantee that the full wave equation, without any approximation, remains 
invariant under the parity-time transformation [Ml [37]. Throughout this work, we 
place our emphasis on the temporal domain and consider a spatially homogeneous yet 
time-dependent dielectric permittivity, e {t). For the sake of illustration we choose 

e (t) = e (r) = d-t-ei cos^ (r)-I-sin (2r) , (1) 

as a modulated time-dependent optical potential [c.f. Fig. [T], similar to its spatial 
counterpart as discussed in Refs. [iniEoiQ. In Eq. ([T|), £ is the background dielectric 
constant, Ei represents the amplitude of the real prohle of the potential, whereas £2 
describes the strength of the gain/loss periodic distribution. Moreover, r = fet is a 
dimensionless time, where b > 0 acts as a scaling factor and indicates the rate (i.e., 
the frequency) of modulation of the permittivity, that we assume to occur slower than 
the oscillations of the wave. This is reminiscent of the similar form of modulation 
adopted in Refs. [Ml [20] for the spatial case and could be experimentally realized by 
utilizing electro-optical systems [3H]- Note that we treat the amplitudes £i and £2 as 
signed quantities in Eq. ([1]), though in general the behaviour is not symmetric under an 
exchange of the sign of the amplitudes. 

f Note that this choice of the frequency-independent permittivity, being the analogue of the refractive 
index used in m in position space, does not satisfy the standard Kramers-Kronig relations, i.e., the 
medium is not subject to the causality principle, e (t) is valid for all time (that is, does not vanish for 
negative times) so that a priori no restrictions on the past are necessary. 
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Figure 1. Symmetric real, 5i(£), and anti-symmetric imaginary, 3(e), parts of the 
time-dependent dielectric permittivity. 


In order to investigate dynamical and spectral properties of light in non-stationary 
PT-symmetric structnres, we derive an exact second order differential eqnation from 
Maxwell’s equations [39] for the electric displacement vector 

X»(r,t) =0, (2) 

which is valid for an arbitrary shape of the time-dependent dielectric permittivity |40] . 
Here, A is the Laplace operator and c is the speed of light in vacuum (for the theory of 
non-stationary electromagnetics see, e.g.. Refs. [m 1121 SSI). Without the modulation 
rate, i.e., when 6 = 0, the standard linear dispersion relation k = 0 J\Js (0)fc/c holds, 
where £ (0) = e + Si is the permittivity at r = 0 and k is the unit vector in the 
direction of propagation. In the presence of modulation ([T]), both the amplitude and the 
angular frequency of light undergo a time-dependent modihcation governed by Eq. (12|) . 
Accounting for this instantaneous effect we seek a solution of Eq. ([2]) by making the 
ansatz 

X>(r,f) = , (3) 

which reflects the spatial homogeneity of the permittivity. Here, it is the unit vector 
along the polarization direction, while the comp lex-valued ‘amplitude’ P describes the 
influence of the modulated potential on the light. In the absence of any modulation, we 
expect to recover the free propagation of light through a uniform medium with a constant 
dielectric permittivity so that P = 1. Note that a similar (full) wave equation for the 
space-dependent electric held and permittivity is discussed in Ref. [36] for describing the 
so-called PT-symmetric coherent-perfect-absorber laser. Moreover, time reversal and 
time-dependent wave propagation is studied in various aspects, such as for time-localized 
perturbations combined with spatial periodicity [S] HS] EG] HT] and sigmoidally changing 
systems with either real or complex permittivity/refractive index [iQl EH] EHl EO] 1^ . 
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Next, we insert the ansatz (E]) in Eq. (|2]) and obtain a second order linear differential 
equation for iF, 


-1 ^(r) =0, 


( 4 ) 


where “dot” refers to the derivative with respect to the dimensionless time r. As our 
interest is restricted to modulations of the complex dielectric permittivity profile, which 
are slow when compared to the oscillations of light, we can adopt fe/cu <C 1 and henceforth 
safely ignore the first term in Eq. (|4]). In this approximation, the remaining first order 
differential equation generally determines the instantaneous angular frequency as 


fl (t) = oj 


(-) 

u 

\ 6(0)3? [£(r)]' 


2 

[ k(r)|^ J 


( 5 ) 


for an arbitrary form of the dielectric permittivity 1§|. The exact solution of the reduced 
equation when integrated from the “initial” time 0 to some instant of time r leads to 


T (r) = exp — ^arctanh(Ci3) 




( 6 ) 


X exp 


uC 


arctanh 


C{B — f^tanr) 


that explicitly exhibits the VT symmetry of the displacement, = T). In Eq. (E]), the 
constant parameters A, B, C are introduced for the sake of brevity: A = e/{e Si) > 0 
carries information about the real potential, whereas B = S 2 /{e + £i) amounts to the 
comp lex-valued permittivity, being the signature of the gain/loss mechanism. Both A 
and combined with C = 1/\fATr~B^ > 0, demonstrate the significance of the real 
and imaginary parts of the permittivity in the instantaneous character of light. Note 
that any modulation vanishes if S = 0, ^ = C = 1 (£i = £2 = 0) and/or b = 0 so 
that we recover the anticipated free propagation of light, as also limb_j.o = 1. The 
solution ([6]) holds for a class of potentials of the form ([1]), which is fully determined on 
choosing A, B and one of the constants in the potential, say e. In addition, to mark 
out the range of variation of these parameters, we insert the (approximate) solution ([6]) 
into the exact Eq. (jll) and estimate the error for a given modulation frequency b and 
angular frequency u numerically. As can be easily checked, the relative error (that is, 
the term {b/u:)‘^ T/T) does not exceed 0.13 for A G [0.7,1.7] and B G [—0.8, 0.8] for a 
ratio of b/u = 0.01 [Jj]. Moreover, other analytical solutions for PT-symmetric (quantum 
mechanical) potentials can be found in Refs. [53l |5l] . 


§ This reminds us of the analogous definition of the local wave vector in the spatial domain (see, e.g., 
Ref. [5^1. 

II The choice for the ratio of modulation and wave frequencies is to estimate the upper value of b/uj, 
for which our approximation is accurate. The effects, as proposed below, remain the same for smaller 
values of b/uj. 
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The solution (Ej) allows us to fully characterize the dynamics of waves in non-stationary 
complex potentials. Indeed, by decoupling the real and imaginary parts of the time- 
dependent component of the electric displacement, (r) = \T> (r)| we obtain 

direct access to the prohle of the instantaneous amplitude \T> {t)\ 

\'^ {'^)f = = exp ci;C&“^arctanh(Ci3)| 


(7) 


X exp 


2B 


uC , 

—-arctanh—--— , 

2b C{A-\- 2B^ tan^ r) 


and also to that of the instantaneous phase 


= 


UT 


uC 
4& V 


arctan 


AC tan r 


arctan ■ 


AC tan r 


1 + BC 1-BC 

the hrst derivative of which yields the prohle ([5]) of the instantaneous angular frequency, 
= —6$, as one would expect. Its explicit form expressed in terms of parameters A 
and B is 


nir) = 


u 


1 + 


A sin^ r -|- cos^ r 


( 8 ) 


(^^sin^ r -|- cos^ -I- B"^ sin^ (2r) 

Again, note that when the VT modulation is “switched off”, the relations \T>\'^ = 1, 
$ = —ut and = a; are obtained (see Ref. [55] for main study methods of signals whose 
frequency content changes in time). 

For a complete description of the instantaneous amplitude and angular frequency, 
and in order to reveal their specihc properties, we determine the extrema of (I7|) and dHj) 
at the stationary points to ^ (“Oo, oo). They are given implictly by the equations 


bA^\ 

2^^ 


\T>f = -Bue^\T>f sin {2 t) = 0, 


T2 = cue^sin (2r) I (1 — [3? (e:)]^ 

(^A sin^ r - cos'^ r) | 

= 0. (9) 

From here we immediately recognize that both quantities have extrema at Tq = 7rm/2 
provided that m is an integer whose even values, m = 2J\f (with J\f being an integer), 
result in |X>(7rA7)|^ = 1 and (ttA/") = oj. In contrast, the odd values, m = 2Af 1, 
give rise to 

|X> (7r/2 -|-ttAA)!^ = exp a;C&“^arctanh(Ci3)| , (10) 

^l{7^/2 + 7^^^) =uj{l + A)/{2A), (11) 

as obtained from Eqs. (I7|) and ([8|), respectively. Equation ffTOj) leads to the maximum 
possible amplitude amplihcation {B < 0) and attenuation {B > 0) correspondingly 
linked to the well-dehned temporal domains of gain and loss. A feature of our complex 
and PT-symmetric potentials is the presence of additional extremal values at 

sin^ r' = (l - SAB^Vd) / (1 - ^) (12) 
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Figure 2. Instantaneous amplitude (left panel, Eq. (O) and angular frequency (right 
panel, Eq. ([5]) or (l8|)) of light modulated by the PT-symmetric dielectric permittivity 
e, Eq. ([T]), for different values of A and B. \B\ < |1 — ^| /2 corresponds to the upper- 
panel (bellow the threshold), \B\ > |1 — ^| /2 to the lower one (above the threshold). 


when the discriminant D = — of the quadratic polynomial (in sin^ r) 

appearing in the curly parentheses in Eq. ([9]) is positive, that is when 2\B\ > |1 — ^|. 
The threshold, however, does not guarantee that the additional extrema are to be found 
at real values of Tq. For this we need to specify the value of B further, albeit the precise 
criterion depends on the value of A. For 0 < ^ < 1, the criterion for additional extrema 
is given by 2|;B| > \/l — A determined by the constraint sin^ Tq > 0, whereas for ^ > 1, 
the criterion changes to 2\B\ > ^A? — A following from the condition sin^Tg < 1. For 
^ = 1, the two criteria are equivalent, and additional extrema can be expected from 
Eq. ([9]) for any real value of B aX t'q = {2M + l)7r/4. 

The physical reason for this difference is the sign of the modulation in the real part 
of the dielelectric permittivity in Eq. ([T]). For 0 < ^ < 1 the modulation of the real 
part is positive, that is, the real part varies between e and e + Si, which gives rise to 
an increase of the instantaneous frequency as shown in Fig. |2b. For ^ > 1 the signed 
amplitude Si is negative and the permittivity varies between e — |£i| and e, which leads 
to a decrease. In both cases the position of the additional extrema is aptly described by 
Eq. ffT2l) . which gives rise to a splitting of the modulation extrema (as shown in Fig. [2li) 
that would otherwise be determined by Eq. flTT]) in the case when the instantaneous 
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angular frequency is modulated by the conventional real permittivity. The value of the 
instantaneous angular frequency at these additional extrema for a modulation by the 
imaginary permittivity beyond the threshold is then given by 

fi(r') = ujC^ (l + 5A + 4:B^ + Vd/ / 8 . (13) 

It is important to note that even though changing the time origin in Eq. ([1]) does affect 
the formal dehnition of time-reversal symmetry, the observed effects remain the same 
with the only modihcation that the symmetry point is shifted for modulations of both 
the instantaneous amplitude, Eq. ([7]), and angular frequency, Eq. (jS]). 

Time-periodic PT-symmetric optical potentials feature unusual, though expected 
modulations of the instantaneous properties of the light. Figure [2] illustrates the 
evolution of these quantities, possessing time-reversal symmetry (|X>(r)|^ = |X>(—r)|^, 
(r) = f2(—r)), against the dimensionless time r for various values of A and B. 
Fig. [2ti shows that the sign of B determines whether the overall dynamics give rise 
to amplihcation {B < 0) or attenuation {B > 0), despite the fact that in both cases 
there are equal periods of gain and loss, occurring, however, with different strengths. In 
addition to these amplitude variations. Fig. [2]3 shows modulations of the instantaneous 
angular frequency towards either higher (0 < ^ < 1) or lower {A > 1) frequencies, 
compared to uj. For small values of \B\, the instantaneous angular frequency remains 
mostly unaltered (Fig. |2 ]d - the additional maxima for ^ = 1 at Tq = (2A/'-|-l)7r/4 are too 
small to discern on this scale) and only the amplitude experiences a modulation. Larger 
values of B lead to a vigorous modulation of the amplitude attenuation (Fig. and 
to a pronounced modihcation of the frequency modulation prohle (Fig. |2ji). Unlike the 
ordinary modulation of angular frequency, where the modulation rate is commensurate 
with the scaling factor (i.e., the modulation frequency b) of the real potential [56], in 
our PT-symmetric structure the extrema of frequency modulation experience a split 
beyond the threshold 2 |S| = |1 — A\, as depicted in Fig. [2li. The global extremum 
turns into a local one and two new global extrema appear on either side such that 
now they occur with a shrunken period and the troughs of curves are shifted towards 
the lower frequency domain, quantitatively determined via Eq. flT^ . The values of the 
instantaneous amplitude and angular frequency at the regular and additional extrema 
are shown in Fig. [3] as functions of A and B, where a comparison is made between 
the VT- and real-potential induced surfaces. As the gain/loss parameter B changes 
its sign from negative to positive, the extrema of the instantaneous amplitude descend 
from the region of amplihcation (|T> (7r/2-|-TrA/")!^ > 1) ^^e region of attenuation 

(|X> (7r/2 -I- 7rJ\f)f < 1) for all values of A (Fig. |3t). By comparison, the unit surface 
indicates the absence of the amplitude modulation when the imaginary part of the 
permittivity is “switched off”, that is for B = 0, which marks the line along which both 
surfaces intersect. In contrast, as seen from Fig. (Sb, the extrema of the instantaneous 
angular frequency differ from unity even if the imaginary part of the permittivity is 
zero. While the S-mdependent surface designated by Q {7i/2 -\- nAf) /u describes the 
extrema of the instantaneous angular frequency as modulated only by the real part of 
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Figure 3. Extrema of the modulus square of the instantaneous amplitude ((a), 
Eq. (fTIill l and angular frequency ((b), Eqs. (ITT]) - (fT^ l. The comparison is made between 
modulations of light via a real permittivity profile (given by the unit surface in (a) 
and the upper surface in (b)) and PT-symmetric permittivity profiles (given by the 
surface crossing unity in (a) and the lower surface in (b)). 


the permittivity, the surface (tq) /oj always lies below (7r/2 + tiM) /oj and represents 


the split extrema due to the imaginary part of e. 

Until now, we have discussed how the VT modulation of the dielectric permittivity 
affects the instantaneous characteristics of light for various values of A and B, as allowed 
for within our approximation. Since the time-dependent permittivity in experimental 
situations can be modulated with the amplitude much less than through the background 
dielectric constant (see, e.g.. Ref. [56]), it is sensible to discuss the impact of the 
permittivity in the dynamics of light for those values of parameters which are currently 
available especially in PT-coupled waveguide devices. If we consider a modulation of 
the dielectric permittivity with the real |£i| <C £ and imaginary \e 2 \ -C d parts and keep 
terms up to hrst order of |£i| /e and \e 2 \ /e m. the parameters A ^ 1 — £i/d, B ^ 62/6 
and C ~ 1 -b £ 1 / (2d), we can reduce Eq. ([7|) to an experimentally accessible form 



(14) 


As seen, the instantaneous amplitude no longer depends on the value of real part of 
the permittivity modulation given by £ 1 , but only on the gain/loss distribution strength 
£ 2 , the sign of which suitably dehnes the temporal domains of gain (£2 < 0) and loss 
(£2 > 0), as also shown in Fig. 0^. As the magnitude of the absolute value of £2 increases, 
the amplitude modulations become more pronounced. For a given |£ 2 | /e = 0,003, for 
instance, one can clearly distinguish between the increase of ~ 35 % and decrease of 
~ 25 % in the amplitude modulation extrema [the black solid and dashed curves]. Such 
an asymmetry in the gain and loss domains is due to the fact that Eq. (1141) changes its 
form under the change of the sign of £ 2 . In contrast to the amplitude modulations, that 
distinctly occur already for small values of £ 2 , the modulations of the lAF are driven 
only by the real prohle since the gain/loss strength contributes with the second-order 
term in the denominator of Eq. ((51). However, if we allow for a strong VT coupling by 
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Figure 4. Instantaneous amplitude ((a), Eq. (fHl) ') and frequency ((b), Eq. (0) of 
light modulated by the PT-symmetric dielectric permittivity e, Eq. o, for the small 
value of the amplitude ei of the real prohle (|ei| /e = 0,001) and for selected values of 
the gain/loss strength £ 2 - 


setting \e2 \ /e < 1, but keep the same restriction for the real amplitude, |ei| /e 1, 
the instantaneous frequency will be modulated beyond the threshold, and therefore, the 
split of extrema and associated shrunken periods can clearly be seen at Tq = (2A/'+l)7r/4 
since A 1, as demonstrated in Fig. Hb- Moreover, when £ 2 /^ = 0.3, the change in the 
instantaneous frequency rises by 4 %, as compared to the initial frequency u. 

4. Conclusions 

In conclusion, we have examined the impact of the time-dependent PT-symmetric 
dielectric permittivity in the dynamical and spectral features of light. In our formalism, 
we have shown that the VT modulation of light is associated with the well-dehned 
temporal domains of gain {B < 0) and loss {B > 0). We have also determined two 
different criteria for the splitting of the extrema of the angular frequency modulation 
to occur and we have demonstrated the shrinkage of the modulation period. Both the 
split and shrinkage are general, inherent features of time-dependent complex potentials. 
A direct manifestation of time reflection symmetry in our particular non-stationary 
structure is always evident. Such effects warrant more detailed future investigations with 
different time-dependent "PT-symmetric potentials. It is true that while the predicted 
curves on Fig. 0] could be observed experimentally with techniques available nowadays, 
the curves on Fig. [2] are likely to only be accessed in future. 

We should mention that our choice of the permittivity suggests that the causality 
principle is not fulhlled, i.e., the parameters of the permittivity do not follow 
conventional dispersion constraints, as represented by the Kramers-Kronig relations. 
Given the known inconsistencies with the Kramers-Kronig relations for PT-symmetric 
systems [571 EH ESI [60] and other artihcial metamaterials with complex-valued 
permittivity [HU ESI |63l[6ll[65], analogous relations must be constructed for frequency- 
independent PT-symmetric time-varying dielectric permittivities, relaxing the strict 
assumptions of the causality. 
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Although the theory developed here can be readily expanded for studying similar 
VT -induced effects for matter [66l [671 EHl IM] and acoustic [5T1 |70l [7T1 [72] waves, it also 
has indirect implications for time-dependent coupling in mechanical systems [731 EU- 
The consideration of the analogous theory for modified VT symmetries [751 [761 [77| 
would be of great interest. 
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